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ABSTRACT 

We report the discovery by M. Linnolt on JD 2455665.7931 (UT 2011 April 
14.29) of the sixth eruption of the recurrent nova T Pyxidis. This discovery was 
made just as the initial fast rise was starting, so with fast notification and re- 
sponse by observers worldwide, the entire initial rise was covered (the first for 
any nova) with fine time resolution and in three filters. The speed of the rise 
peaked at 9 mag/day, while the light curve is well fit over only the first two days 
by a model with a uniformly expanding sphere. We also report the discovery by 
R. Stubbings of a pre-eruption rise starting 11 days before the eruption, peak- 
ing 1.1 mag brighter than its long-time average, and then fading back towards 
quiescence five days before the eruption. This unique and mysterious behavior 
is only the fourth known (with V1500 Cyg, V533 Her, and T CrB) anticipatory 
rise/dip closely spaced before a nova eruption. We present 19 timings of pho- 
tometric minima from 1986 to February 2011, where the orbital period is fast 
increasing with P/P = -^313, 000 years. From 2008-2011, T Pyx had a small 
change in this rate of increase, so that the orbital period at the time of erup- 
tion was 0.07622916±0. 00000008 days. This strong and steady increase of the 
orbital period can only come from mass transfer, for which we calculate a rate of 
-|^g-6.o±o.5 yT~^ ioT mass leaving the companion star. We report 6114 magni- 
tudes between 1890 and 2011, for an average B = 15.59 ± 0.01 from 1967-2011, 
which allows for an eruption in 2011 if the blue flux is nearly proportional to the 
accretion rate. We present a model for the infrared-to-X-ray emission during qui- 
escence as being from a 34,000 K blackbody nearly filling the Roche lobe (caused 
by the very high accretion rate creating an extended envelope around the accre- 
tion column) plus a z/°'^ non-thermal component (from the optically-thin outer 
and circumbinary regions). 

Subject headings: novae, cataclysmic variables - stars: individual (T Pyx) 
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1. Introduction 

T Pyxidis (T Pyx) was the first discovered recurrent nova (RN) with three known 
eruptions, and this over a dozen years before any other RN had even their second discovered 
eruption (Schaefer 2010). T Pyx is one of the quintessential RN, with six eruptions in 1890, 
1902, 1920, 1944, 1967, and now 2011. Over the past few decades, the quiescent magnitude 
of T Pyx has been near V=15.5, while peaking at V=6.4 (Schaefer 2010). T Pyx is slow 
among RNe, with its overall rise taking 40 days and the decline by three magnitudes from 
the peak {t^) taking 62 days (Schaefer 2010). T Pyx is one of two RNe with an orbital 
period (P=0.076 days; Schaefer et al. 1992; Patterson et al. 1998; Uthas et al. 2011) that 
is inside or below the period gap. T Pyx is the only RN that has a nova shell, and this 
shell has a surprisingly slow expansion velocity plus a structure of thousands of discrete 
knots (Duerbeck & Seitter 1979; Wilhams 1982; Shara et al. 1989; 1997). The homologous 
expansion of the individual shell fragments proves that they suffer no significant deceleration, 
and the expansion velocities, viewed over 13 years by Hubble Space Telescope, prove that the 
nova shell was ejected in the year 1866±5 at a velocity of 500-715 km s~^ (Schaefer et al. 
2010a). This velocity is greatly too low to be from a RN event, and with the measured mass 
in the shell (~10~^'^ M0) being much greater than is possible for a RN event, we know that 
the 1866 event must have been a normal classical nova eruption with a very long time of 
prior quiescence (Schaefer et al. 2010a). 

A key and unique fact about T Pyx is that its quiescent brightness has been secularly 
fading from 1890 to 2011 (Schaefer 2005). Before the 1890 eruption, archival plates show T 
Pyx to be constant at B=13.8, the next four inter-eruption intervals have B-band magnitudes 
of 14.38, 14.74, 14.88, and 14.72, while after the 1967 eruption the quiescent B magnitude 
fades from 15.3 (around 1968) to 15.5 and fainter (from the 1980's until 2011). This secular 
decline is fundamental for understanding the evolution and fate of T Pyx, as well as the 
timing of the eruptions. 

T Pyx has long been a perplexing system, because there was no known way to combine 
a short orbital period and a high accretion rate into one system. (The accretion rate must 
be high, ~10~^ Mq yr~^, so as to provide enough mass to trigger eruptions with a typical 
recurrence time interval of two decades.) That is, a binary with a period below the gap 
should have angular momentum loss only from gravitational radiation, and this will lead to 
an accretion rate many orders of magnitude smaller than required (Patterson 1984). This 
dilemma was solved by Knigge et al. (2000) who realized that the high accretion rate could 
be powered by a hot source on the white dwarf that is irradiating the companion star. For 
T Pyx, this idea was extended by Schaefer et al. (2010a), where the classic nova eruption 
of 1866 started a high accretion rate which drives a largely self-sustained supersoft source 
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(with nuclear burning on the surface of the white dwarf). This mechanism can only work 
for binaries with short orbital periods and where a substantial magnetic field channels the 
material into a small area. The secular fading demonstrates that the supersoft source is not 
fully self-sustaining. Schaefer & Collazzi (2010) generalizes this situation to a whole class 
of stars, called V1500 Cyg stars, which suffer this same secular fading after their eruptions, 
and which also have short orbital periods and highly magnetized white dwarf stars. So 
the picture we have for T Pyx is that it started out as an ordinary cataclysmic variable 
with a long period of quiescence, then an ordinary nova eruption in 1866, which ignited a 
supersoft X-ray source whose irradiation drives a high accretion rate and hence recurrent 
nova eruptions. Over time, the accretion rate declines and the quiescent brightness fades. 

By far the most important question relating to RNe is whether they are the progenitors 
of Type la supernovae. RNe are one of the best progenitor candidate classes, because they 
certainly have white dwarfs near the Chandrasekhar mass being fed with a very high accretion 
rate. Nevertheless, T Pyx (and any short-period RNe) cannot be a Type la progenitor 
(Schaefer et al. 2010a; Schaefer 2010). The reason is that the ordinary nova eruption in 
1866 easily dominates the dynamics of the system, and ordinary nova eruptions eject much 
more mass than has been accreted by the white dwarf, so the white dwarf is actually loosing 
mass and will not become a supernova. Another way of saying this is that the secular 
fading shows that the interval during which T Pyx is a RN is short, perhaps under two 
centuries, and thus will have httle effect even if the RN event causes a net gain of mass on 
the white dwarf. The long-period RNe remain as strong progenitor candidates (due to their 
high accretion rate being driven by the evolution of the donor star so as to be sustained and 
long in duration), but the short-period RNe cannot sustain the high accretion and so are 
not supernova progenitors. 

The secular fading is also critical for understanding the timing and triggering of the 
RN events. Schaefer (2005) demonstrated that the inter-eruption intervals during which 
the quiescent magnitude was faint on average lead to a long time between eruptions, while 
a bright interval lead to a short time between eruptions. The physics is simple, with the 
amount of accreted mass required to trigger an eruption being constant for any one RN, with 
the accretion rate being given by the blue flux dominated by the accretion disk, so that a 
high quiescent brightness implies a high accretion rate and a short interval needed to accrete 
the trigger mass. With calibration from prior inter-eruption intervals, it becomes possible 
to predict the approximate date of upcoming eruptions. As such, Schaefer (2005) predicted 
that U Sco would erupt in 2009.3±1.0, and this advance notice allowed for an incredible 
coverage of the eruption in 2010.1, so that this now is the all-time best observed nova event 
(Schaefer et al. 2010b). Schaefer (2005) also predicted that the next eruption of T Pyx 
would be in 2052±3, while Schaefer et al. (2010a) extended this to predict that the next 
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eruption will be after the year 2225. This prediction was shown to be horribly wrong by the 
T Pyx eruption in 2011, so the issue now is to understand where the prediction went wrong. 

2. Discovery of the 2011 Eruption 

Starting with the 1890 eruption, the inter-eruption intervals were 11.9, 17.9, 24.6, and 
22.1 years (Schaefer 2010). In the middle 1980's many people around the world (including 
many of the authors of this paper) realized that T Pyx 'should' go off any year now. The 
simple logic was that the prior inter-eruption interval was 22 years, and 22 years after the 
1967 eruption placed the expected eruption around 1989. Due to this, many of the world's 
best observers kept T Pyx under nightly watch. 

As the 1990's stretched on, the expected eruption did not occur. Vigilance was main- 
tained and occasional public exhortations were published in many journals. Any eruption 
after 1967 would certainly have been discovered because the T Pyx eruption duration is two- 
thirds of a year and so no eruption would have been missed due to observational gaps. By 
the year 2005, the interval since the prior eruption extended to 38 years, with no explanation 
for this long delay. 

Schaefer (2005) provided the explanation. The key point was that T Pyx has been 
declining overall in brightness ever since before the 1890 eruption, and this means that 
the accretion rate in the system was declining, so that it would take longer to accumulate 
the trigger mass on the surface of the white dwarf. In particular, from the 1944 to the 1967 
eruption the average B-band magnitude was 14.72±0.03, while from 1967 to 2005 the average 
B-band magnitude was 15.49±0.01 mag. The blue flux from T Pyx is entirely dominated by 
the accretion disk, so the B-band flux is determined by the accretion rate. The dimming by 
0.77 mag (a factor of two drop in accretion luminosity) requires a fall in the accretion rate 
by at least a factor of two, and this forces an increase in the inter-eruption time interval. 
The long delay since 1967 was simply because the rate of accumulation on the white dwarf 
surface had at least halved. 

With the prediction from Schaefer (2005) that the next eruption was a long way in 
the future, professional preparations for the next eruption (like the long-running Target of 
Opportunity program with the Hubble Space Telescope) were stopped. But this did not deter 
many amateur observers (who were generally fully aware of the prediction) from continuing 
to regularly check for eruptions, report visual magnitudes, and make long photometric time 
series with CCDs. For example, Stubbings started in early 1998 in making 2004 visual 
estimates of the brightness of T Pyx with a 16-inch reflecting telescope in Victoria Austraha, 



-6- 



while Richards made CCD time series observations on five nights in early 2008 with 661 
brightness measures. 

Prom 2008 until 31 March 2011 (JD 2455652), T Pyx held at a steady magnitude of 
V=15.5, with the usual periodic variations of <0.2 mag. Por the 610 observations in the last 
4000 days of this interval, Stubbings never saw T Pyx brighter than V=15.0. Thus, it was a 
distinct surprise when Stubbings saw T Pyx at V=14.5 on JD2455657.0104 (UT 2011 April 
5.51), with the realization that this could be an initial rise of the next T Pyx eruption. With 
this exciting possibility, Stubbings checked T Pyx five more times before dawn, but saw only 
variations from 14.4 to 14.7. On the next night T Pyx was slightly fainter, and it continued a 
slow fade until JD2455662.0417 (UT 2011 April 10.54). The light curve just before the sharp 
rise to the eruption (see Figure 1) shows a significant short-duration pre-eruption rise, which 
started roughly 13 days before the eruption, peaked at around V=14.4 mag (close to 1.1 
mag above the prior quiescent level), and declined towards quiescence before the eruption. 

The existence of this pre-eruption rise has been confirmed with the archival data from the 
All Sky Automated Survey telescope located at the Las Campanas Observatory in northern 
Chile (Pojmanski 2002). The images were taken through a Bessel V filter and magnitudes 
are tied to a large number of comparison stars in the field with magnitudes in the Tycho 
catalog. This procedure can lead to moderate offsets for an object with unusual color (like T 
Pyx), and indeed the average magnitude in the year before eruption (15.12) is brighter than 
the average of other V-band measures (V=15.51) over the same time interval. Nevertheless, 
the differential photometry is excellent to within the statistical uncertainties (typically 0.09 
mag). For the year prior to JD 2455654.5972, T Pyx was seen to vary about its average with 
an RMS scatter of 0.14 mag. On JD 2455656.62097 (0.4 days before Stubbings discovery 
of the rise), T Pyx was 0.43 mag brighter than the average, with this being the brightest 
magnitude by far in the prior year. This is the confirmation of the pre-eruption rise. On 
JD 2455658.62564 and 2455660.62294, T Pyx was consistent with being back to the prior 
average. Thus, the apparent duration of the pre-eruption rise is roughly four days with the 
beginning around 11 days before the start of the eruption. 

This pre-eruption event is one of an unusual phenomena that is still completely myste- 
rious. CoUazzi et al. (2009) report an extensive survey of archival data for 22 novae with 
well-sampled pre-eruption light curves, which showed that two of these (V1500 Cyg and V533 
Her) have pre-eruption rises with different time scales and morphology, while one nova (T 
CrB) has a significant pre-eruption dip. Table 1 summarizes the properties of these events. 
In all cases, the pre-eruption rise cannot be some early part of the usual thermonuclear run- 
away because the anticipation times and the durations (weeks and years) are greatly longer 
than both the theoretical runaway timescale (hours) and the observed rise times (hours). 
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In all cases, the pre-eruption event is highly significant and unprecedented for long time 
intervals before the eruptions and after the eruptions. For example, the dip before the 1946 
eruption of T CrB was deepest (by roughly 2 magnitudes below the usual quiescent level) 
just 29 days before the eruption, whereas there are no other dips from 1867 to 1945 and from 
1947 to 2011. The uniqueness and the time coincidence immediately before the eruption pro- 
vide proof of a causal connection. It is difficult to understand this connection, because the 
pre-eruption brightness level is dominated by the accretion which is governed by the mass 
falling off the donor star, while the timing of the eruption is governed by conditions at the 
bottom of the accretion layer deep under the surface of the white dwarf. How could the 
donor star near the LI point possibly anticipate that conditions deep in the white dwarf will 
soon heat to the trigger threshold? We can only speculatively invent several possibilities. 
Perhaps some isotope within the accreted layer starts burning so as to heat the bottom of the 
layer up to the point where the thermonuclear runaway can start, but we do not see how the 
very small heated surface of the white dwarf can significantly brighten the whole system, nor 
do we know any suitable isotope. Perhaps some random increase in the accretion rate will 
heat up the accretion layer sufficiently so as to trigger the thermonuclear runaway, but the 
pre-eruption increases are so unique (especially that of V1500 Cyg) that it seems they are 
not from any independent original cause, the extra energy of accretion is too small to make 
any significant rise in temperature at the base of the layer, and the accretion energy cannot 
get to the bottom of the layer on the time scale of days as required for T Pyx. In all, we can 
think of no plausible explanation for the pre-eruption rises, while the pre-eruption dip for 
T CrB does not even have these poor speculations to explain why its turn-off of accretion 
anticipated its eruption. We regard the explanation of these pre-eruption phenomena as a 
premier challenge for nova theory. 

Linnolt had also been monitoring T Pyx since March 2009 from Hawaii, waiting for 
an eruption, as part of a larger program of monitoring recurrent novae for eruptions (e.g., 
Schaefer et al. 2010b). He used a 20-inch aperture Newtonian reflector telescope operat- 
ing visually at 500x magnification, so that observations down to V=17 were possible. On 
JD2455664.8590 (UT 2011 April 13.36), Linnolt observed T Pyx at V=14.5 and realized that 
this was unusually bright, and hence a possible indication of an eruption. With anticipation, 
Linnolt checked T Pyx early the next night (JD2455665.7931, UT 2011 April 14.29), saw 
that the star was bright (V=13.0), and immediately realized that an eruption had started. 
Within 15 minutes, the discovery was announced on the AAVSO discussion group which 
is closely monitored by interested amateurs and professionals worldwide. The discovery 
was quickly confirmed by Australian visual observers Plummer (JD2455665.8847 V=12.2) 
and Kerr (JD2455665.941 V=11.3). CCD observers Nelson, Carstens, and Streamer started 
taking time series photometry in multiple bands by JD2455665.9973. 
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To the best of our knowledge (e.g., Payne-Gaposchkin 1964; Robinson 1975; CoUazzi 
et al. 2009; Strope, Schaefer, & Henden 2010; Schaefer 2010), the 2011 eruption of T Pyx 
has by far the best observed initial rise from quiescence. Only a handful of old novae have 
had just a few points taken during the later half of their initial fast rise, whereas T Pyx has 
high time resolution time series with near full-time coverage and in multiple bands. Indeed, 
T Pyx is the only nova with any observations in the middle or the first part of the initial 
fast rise. The very fast alert to the world and the very fast reaction by observers worldwide 
made this possible. 



3. The 44 Year Inter-Eruption Interval 

Schaefer (2005) and Schaefer ct al. (2010a) made a quantitative physics-based prediction 
that the date of the next T Pyx eruption would be in the year 2052 or much later. The first 
key idea is that the accretion rate (and hence the accumulation rate on the white dwarf) 
can be measured from the blue flux, so a B-band light curve between eruptions can be used 
to determine the average accretion rate. The general relation will be 

M = MieFf = MielO-"-^^^^-^'^) , (1) 

where M is the accretion rate, Miq is the accretion rate when T Pyx has a B magnitude of 
16.0, B is the B-band magnitude, Fb is the B-band flux in units of the the flux when B=16, 
and E an exponent expressing the physics of the light sources in T Pyx. The use of the 
B=16 mag level is arbitrary (any other magnitude level could have been adopted), but this 
choice is convenient for T Pyx. The second key idea is that the total mass accreted between 
eruptions can be estimated from the accretion rate averaged over the inter-eruption interval 
of AT, with 

AMacc = {M)AT. (2) 

Here, the accretion rate values from equation 1 are averaged over all available blue magni- 
tudes. The third key idea is that the trigger mass will be close to the mass accreted between 
eruptions: 

Mtrigger = AMacc- (3) 

The fourth key idea is that a given recurrent nova system should have a unique trigger mass 
that will be the same for all eruptions. The trigger mass can only be a function of the white 
dwarf mass, radius, magnetic field, and the composition of the accreted material, and these 
vary from RN to RN, but not from eruption to eruption. This can be expressed as 

^triggerl — ■^trigger2 — ^triggerZ — • • • i (4) 
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where the numerical subscript indicates different eruptions. Putting these equations together, 
we see that (Fj')AT should be a constant from eruption to eruption. For intervals where 
the RN is bright on average, then the accretion rate will be high, the trigger mass will 
be accumulated relatively fast, and the inter-eruption interval will be relatively short. For 
intervals where the RN is faint on average, the time between eruptions will be long. For 
both T Pyx and U Sco (Schaefer 2005; Schaefer et al. 2010b), the blue flux was seen to vary 
with the short inter-eruption intervals having a bright quiescence (and thus a high accretion 
rate) and with the long inter-eruption intervals having a faint quiescence (and thus a low 
accretion rate). We now can determine the current inter-eruption time interval as 

AT3 = AT2(10-°-^^(^-^6))2/(10-°-^^(-^-^6))3. (5) 

The subscripts on the brackets indicate an average over the indicated inter-eruption interval. 

We now have a means of predicting the date of the next eruption by simply measuring 
the historical light curve along with some assumption for the continuing brightness. By 
calibrating with previous inter-eruption intervals, the various constants of proportionality 
cancel out, and we need never know the trigger mass, accretion rate, or the distance. Schaefer 
(2005) predicted that the next eruption of U Sco would be in 2009.3±1.0, and this formed 
the basis for a large international collaboration of researchers at all wavelengths designed to 
discover and then intensively observe the eruption. With the discovery of the eruption in 
2010.1 by B. Harris and S. Dvorak, this advance preparation allowed the U Sco eruption to 
become the all-time best observed nova event. Schaefer (2005) also predicted that T Pyx 
would next erupt in 2052±3, while Schaefer et al. (2010a) extended this date to the year 
2225 or much later. 

This prediction for T Pyx was refuted by the 2011 eruption. So the task now is to 
understand why the prediction went so far wrong. The basic logic of the prediction is simple 
and easy. The light curve is confident and the mass for the trigger threshold really has 
to be constant. About the only place that could be wrong is the conversion from the B- 
band magnitude to an accretion rate. Schaefer (2005) calculated that the accretion rate is 
proportional to the blue flux raised to the 2.0 power (i.e., E = 2), with this exponent being 
derived from a model for a simple accretion disk with the parameters for T Pyx. (For U Sco 
with its much larger disk, the exponent is 1.5.) So this exponent is the likely cause of the 
wrong prediction. 

If the exponent is around 1.0, then the predicted date would have been reasonably close. 
This is the case where the accretion rate is simply proportional to the B-band flux. The 
quiescent magnitude from 1944 to 1967 was 14.72, while from 1967 to 2005 it was 15.49 
(Schaefer 2005), for a magnitude difference of 0.77 mag, corresponding to a ratio of B-band 
flux of close to 2.0. For the accretion rate being proportional to the B-band flux, the average 
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accretion rate must then also have a ratio of close to a factor of 2.0, and the post-1967 interval 
must be 2.0 times that of the pre-1967 interval. This would imply an eruption 2.0x22 years 
after 1967, which predicts the next eruption in 2011. So it seems that, somehow, T Pyx 
effectively has the B-band flux being proportional to the accretion rate (i.e., E 1). 

This fast and easy analysis can be substantially improved. Using the quiescent light 
curve and the inter-eruption time intervals, the value of E that makes for the most constant 
(FJ')AT from interval to interval can be found. 

To this end, we have accumulated a large number of magnitudes of T Pyx in quiescence. 
Table 2 presents these magnitudes, although the collected limits and visual estimates are 
not included. (The printed table only shows only the first 10 lines for illustration, while 
the full table is available in the electronic paper.) The total is 6114 magnitudes, including 
3856 magnitudes during quiescence from 1890-2011, 15 magnitudes during the pre-eruption 
rise just before the 2011 eruption, and 2243 magnitudes in the first two days of the 2011 
eruption. Figure 58 of Schaefer (2010) already displays most of the B-band magnitudes, 
showing the secular dechne since 1890. The magnitudes before 1954 are all from the Harvard 
College Observatory plate archives and are tabulated here for the first time. After the 1967 
eruption, magnitudes were reported by Mumford (1971), Landolt (1970; 1977), Schaefer et 
al. (1992), and Patterson et al. (1998). We also report many more magnitudes for the 
first time, including many individual nights from 1987-2009 with BVRI photometry with 
various telescopes at Cerro Tololo Inter- American Observatory by Schaefer, time series on 
five nights in 2008 at Pretty Hill Observatory by Richards, time series on three nights in 
2011 in Huntington California by Myers, and time series on five nights in 2011 at EUinbank 
Observatory by Nelson. The time series in 2011 were made in support of observations by 
the Chandra X-ray observatory. By good fortune, these time series were as close as 48 days 
before the eruption, and this is critical for getting an accurate pre-eruption orbital period. 
We have collected more unpublished magnitudes that are recorded in the AAVSO database. 
Many more magnitudes in quiescence are already known (mainly in the large database of the 
Center for Backyard Astronomy as reported in Patterson et al. 1998 and Uthas et al. 2011) 
and can be measured (mainly with the archival plates from the 1920's to the 1980's at the 
Sonneberg Observatory), with these extra magnitudes able to fill some of the gaps in time. 

To determine the best fit exponent, we have taken the B and V magnitudes from each of 
the intervals. (The R and I magnitudes are few in number and usually nearly simultaneous 
with a B and V magnitude.) The V magnitudes are converted to B magnitudes by using 
the long term average oi B — V = -1-0.12 (Schaefer 2010). Nightly averages are formed, so 
that a few nights with many magnitudes as part of a time series (which might or might 
not be representative of the whole interval) do not dominate the average. The B-band 
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magnitudes are then converted to flux with Equation 1, where E is held as a free parameter. 
These nightly averages are then straight averaged for each inter-eruption time interval. The 
uncertainty is taken to be the RMS scatter divided by the square root of the number of 
nights. (If the nightly averages are first formed into yearly averages and these are averaged, 
then we get fluxes which are essentially identical.) Table 3 gives the pre-eruption intervals, 
their durations (AT from Schaefer 2010), the average B-band magnitudes over each interval 
((5)), the average Fg for those intervals, the values of (Fb)AT (which should be nearly 
constant), and the values of (F^-^^)AT. 

The only free parameter in this analysis is E. The best estimate for E will be that value 
which returns the most constant values for (Fj')AT across all measured time intervals. The 
values and the scatter will scale as E changes, so we find the value of E that has the least 
fractional scatter. The smallest fractional scatter is near E — 1.13, with 13% RMS variations 
across the five intervals with measures. With this. Table 3 has the value of {F^-^^)AT in the 
last column, with a straight average of 71. Table 3 also lists the values of (Fb)AT (with a 
straight average of 64) because later we will provide a theoretical justification that E — 1. 

The scatter in {F^)AT is much larger than the formal uncertainties, and this indicates 
that our basic model is not exact. So, at the 13% level, the imperfect sampling of the hght 
curve might be misrcpresentative, the value of E varies, or the trigger mass might vary from 
eruption to eruption. For whatever reason, this means that the total uncertainty in (Fj')AT 
is around 13%, so for five independent measures the uncertainty on the average is around 
6%. For predicting any single intcr-cruption interval, the uncertainty will be 13%. Adopting 
a 13% uncertainty, we can make a chi-square fit for E, with this giving a value of 1.13±0.14. 
As such, we see that the five inter-eruption time intervals are consistent with E = 1. 

The observed intcr-cruption intervals and light curves are consistent with E = 1.13 ± 
0.14, but we still need to explain the low E value. For this, the basic calculation of Schaefer 
(2005) was repeated with various physically motivated differences. The basic model adopted 
the usual thin a-disk, extending from near the surface of the white dwarf to 90% of the Roche 
lobe radius, with the ordinary disk light (expressed in Equations 5.45 and 5.49 of Prank et al. 
2002) dominating the system. New shows that T Pyx is an intermediate polar and that this 
is driving bright continuing nuclear burning on the surface of the white dwarf (Schaefer & 
Collazzi 2010). We know that T Pyx is an intermediate polar because it displays the typical 
unstable photometric periodicities near the orbital period (Patterson et al. 1998), with this 
being the observational hallmark of intermediate polars. This new evidence means that the 
disk is truncated by the magnetic field, with an inner edge far from the white dwarf, and 
the disk might be only a small annulus or there might not even be any recognizable disk. 
This also means that T Pyx has a bright supersoft X-ray source on the white dwarf that is 
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irradiating the disk, accretion column, and companion star, hence contributing substantial 
luminosity that is reprocessed into the B-band. This supersoft source hs faded over the 
last century, with its luminosity falling as the accretion rate falls. So we now have two 
modifications to the basic model (truncated disks and reprocessing of the supersoft light) 
that might change the theoretical value of E. 

To evaluate the exponent E in M oc for the relevant conditions, hte system brightness 
as a function of the accretion is calculated. The luminosity of the supersoft source will be 
O.OOTXMc^ for hydrogen burning, with X being the hydrogen fraction taken to be the solar 
value of 0.7. The temperature of the reprocessed light is given by Equation 5.93 of Frank 
et al. (2002). The blue flux from the disk comes from the sum of the blackbody brightness 
of each annular ring throughout the disk with the temperature varying with radius. With 
this, we confirm that £■ fa 2 for the case with no reprocessing and no truncated disk. This 
exponent can be varied somewhat with the accretion rate, but even with a low accretion 
rate of 1 x 10"^ yr"^ E is still 1.82. Further, we find the same result with varying 
inner disk radius and reprocessing sites. That is, ~ 2 is still the best model value. The 
lowest exponents come from the inclusion of reprocessing off the companion star, the largest 
radii of the inner edge of the disk, and lowest accretion rates, but we never get E < 1.5 for 
any plausible conditions. A further problem with these extreme conditions is that the total 
luminosity from the system becomes much fainter than observed, even if the distance to T 
Pyx is as close as possible. Under no plausible conditions can we get E approaching unity. 
That is, our more comphcated and reahstic models cannot account for the E — 1.13 ± 0.14 
value required to explain an eruption in 2011. 

We can understand this general result with a simple argument. First we have the lu- 
minosity of all the various components in the system being proportional to the accretion 
rate (L oc M). This includes the luminosity from viscous heating in the disk, the accre- 
tion column, the supersoft source, and all the reprocessed light that is proportional to the 
white dwarf luminosity. (This is not exact, as the disk and the reprocessed light are also 
proportional to the disk height which scales as M*^'^^. Importantly, all the light sources in 
the system that are nearly constant, like the unheated companion star and the white dwarf 
itself, are negligibly faint.) Second, all of the significant light sources in the system have 
luminosity coming from blackbody emission, or as the sum of blackbodies with a range of 
temperatures, for which L oc crT^. Third, for the case of relatively hot blackbodies, we are 
near the Rayleigh- Jeans limit for which Fb oc T. Putting these three scahng laws together, 
we get M (X F^, which is to say E — A. Even with the scahng with disk height, we only 
get E = 4/1.15 = 3.5. In reality, the emitting objects are not so hot as to have the blue 
light being far out on the Rayleigh- Jeans tail, so Fb becomes proportional to a higher power 
of T and E lowers to around 2. The temperatures predicted in the T Pyx system are all 
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relatively high because the system has a very high accretion rate, the supersoft source is pow- 
ered by nuclear burning, and the system is small as indicated by the short orbital period. 
(In contrast, most other cataclysmic variables have much cooler emission from the accreting 
material because they have greatly lower accretion rates, a white dwarf whose luminosity is 
powered only by accretion, and a larger disk with cooler outer parts.) So we have a general 
conclusion that blackbody emission from all the bright system components have E 2 and 
certainly nowhere near £^ ~ 1. 

This general conclusion does not apply if the dominant light source does not emit as a 
blackbody. At first look, a good reason to take the system light to be from a thermal disk is 
that its spectral energy distribution is a rising power law reminiscent of the F^, oc u^^^ that 
is characteristic of accretion disks (Schaefer 2010). But this simple argument fails because 
the power law stretches from K-band to U-band, and this requires a much broader thermal 
disk than is possible for T Pyx. To get a v^/^ spectral region stretching over almost an order 
of magnitude in frequency from a thermal disk requires that the radius of the outer edge of 
the disk be more than a thousand times the radius of the inner edge (see Figure 5.2 of Frank 
et al. 2002). This cannot be the case for T Pyx because the outer edge is truncated by the 
Roche lobe (10^°'^ cm) while the inner edge is truncated by the white dwarf (10^'^ cm) and 
the magnetic field (10^ cm to nearly lO^'''^ cm). On looking closer, the rising power law is 
Fj, oc v^-^ (Gilmozzi & Selvelli 2007; Schaefer 2010), which is substantially bluer than any 
known thermal disk in a cataclysmic variable. (A small portion of a blackbody spectrum 
can have such a slope, but not over the large range from the infrared to ultraviolet.) This is 
emphasized by the incredibly blue U-B color of -0.95 (Mumford 1971; Landolt 1970; 1977) 
or -1.25 after reddening correction (Webbink et al. 1987). (Indeed, we have just recognized 
that the U-V color for T Pyx had been incorrectly transferred from Table 25 to Table 
32 of Schaefer 2010 without the negative sign, and so the spectral energy distribution in 
Figures 70 and 71 show a false turn down in the ultraviolet. Actually, when done correctly, 
the placement of the U-band point is in line with the rest of the power law shape.) The 
overall shape of the T Pyx spectral energy distribution has no precedents among cataclysmic 
variables or related objects (Dubus et al. 2004; Harrison & Gehrz 1992). Thus, the light 
from T Pyx is not from a thermal disk (with or without any of the effects like reprocessing) 
because its spectral energy distribution is unlike any other cataclysmic variable disk and 
because the observed power law has a range in wavelength that implies a very large span of 
disk radii that is impossible for T Pyx. 

The spectral shape is not possible for any blackbody (or sum of blackbodies) within the 
Roche lobe of the white dwarf. The problem is that the spectrum shows no sign of turning 
over to a Rayleigh- Jeans slope {Fy oc z/^) even far out into the infrared, and this makes it 
impossible to produce the optical and infrared flux from any combination of blackbodies 
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within the Roche lobe. The calculated K-band brightness for a uniform blackbody with the 

radius of the Roche lobe and a temperature of 1350 K (to place the spectral peak in the 
K-band) is 23.8 mag, as opposed to the observed K=14.18 (Schaefer 2010). We can match 
the observed brightness only with an emitting surface area roughly 7000 times larger than 
the Roche lobe. If we try adding in a hotter blackbody, then we should start to sec the 
Rayleigh- Jeans tail. For a 3000 K blackbody (to have the Wien peak at 965 nm so that 
the Rayleigh- Jeans tail should be readily visible in the K-band) from the surface area of 
the Roche lobe, we get K=20.8, which is still several hundred times too faint to explain the 
observed flux. If we add in a 34,000 K blackbody with the size of the Roche lobe, then we 
still get only K=16.9. (An effective temperature of 34,000 K has been suggested to account 
for the observed ultraviolet spectrum by Gilmozzi & SelvcUi 2007.) The addition of hot 
blackbodics then results in a mismatch with observations because the required Rayleigh- 
Jeans tail is not seen. No possible combination of blackbodics within the Roche lobe can 
reproduce the observed infrared flux. Even in the R-band, the observed flux cannot be 
accounted for by a blackbody any cooler than 34,000 K, and this is rejected by the utter lack 
of a Rayleigh- Jeans tail in the spectral energy distribution. The B-band flux can be produced 
by a 34,000 K blackbody from a region about half the size of the Roche lobe, but again this 
possibility is rejected by the lack of a Rayleigh- Jeans tail. The addition of blackbody flux 
from a hotter region would not solve the problem, while the addition of blackbody flux from 
a colder region could never add signiflcant flux to hide the Rayleigh- Jeans tail. In all, the 
observed spectral energy distribution from the K-band to the B-band cannot arise from a 
blackbody (or sum of blackbodics) from within the Roche lobe. 

So where does the optical and infrared light come from? One possibility is that the 
light comes from a circumbinary disk. Such a disk could form from material ejected during 
the accretion process, and this material could be quite massive. Any such disk must reside 
substantially further away from the white dwarf than the companion star, and hence its 
material must be relatively cold. This could easily account for the infrared flux. Any such 
disk should be blown away during each recurrent nova eruption, so a test of this circumbinary 
disk possibility will be to look for the infrared flux after the 2011 eruption has completely 
stopped. However, for any plausible temperature, the circumbinary disk cannot account 
for the B-band flux. The only other possibility is that the optical light is from nonthermal 
emission, which is to say from neither optically thick sources nor blackbody sources. Possible 
mechanisms for the optical continuum include cyclotron, bremsstrahlung, and Thompson 
scattering of some interior light. So we conclude that the optical continuum from T Pyx 
must be dominated by some nonthermal emission mechanism, while the infrared light might 
be dominated by some circumbinary disk. 

If the optical light cannot be dominated by blackbodics, then the earlier models (with 
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E 2) all become irrelevant. There are three reasons why the usual a-disk models won't 
work for T Pyx. First, T Pyx is an intermediate polar so high magnetic fields will destroy 
much of the thermal disk structure and allow for cyclotron radiation. Second, thermonuclear 
burning might still be occurring on the surface of the white dwarf, and this emission will 
dominate any disk, with this being rare among cataclysmic variables. Third, T Pyx has 
a very high accretion rate, as it must have to be recurrent, with the local accretion rate 
on the pole caps of the white dwarf leading to unusual configurations (Nomoto 1982), with 
this being rare among cataclysmic variables. T Pyx exhibits all three conditions and thus 
comparison with cataclysmic variables has limited relevance. T Pyx has high luminosities 
and magnetic fields, conditions which favor nonthermal emission. Indeed, the spectral energy 
distribution shows that the optical light is dominated by non-thermal, rather than thermal, 
emissions. 

Nonthermal emission must be optically thin. In various plausible situations, the emitted 
fiux will be proportional to the accretion rate. For example, for cyclotron radiation in an 
optically thin region, the flux will be proportional to the number of electrons, which will be 
proportional to the accretion rate, so that E=l. Or if the geometric volume of the emission 
region scales as the accretion rate (for example, a doubling of the rate leading to a doubling 
of the volume), then even emission by bremsstrahlung (which scales as the square of the 
electron density) will yield E=l. Or if the thermal energy of the electrons in the optically 
thin region is controlled by the luminosity of some underlying optically thick region that 
is proportional to the accretion rate, then the fiux from cyclotron is proportional to the 
electron thermal energy and hence to the accretion rate, so that E=l. 

So we now have a plausible explanation for why T Pyx erupted in 2011. The observed 
fiux from T Pyx has been dechning for the last 121 years, so the accretion rate has been 
falling and it takes longer times between eruptions to accumulate the nova trigger mass. In 
the 22 year inter-eruption interval before the 1967 event, T Pyx was nearly a factor of two 
times brighter on average compared to the next inter-eruption interval, so the next inter- 
eruption interval should be roughly twice as long (i.e., around 44 years). Therefore, the 
next eruption after 1967 should occur in 1967+44, which is to say in the year 2011. A more 
accurate postdiction is to use the average value of (Fb)AT = 64±9 and (Fg) = 1.47 ±0.011 
to get AT = 43.6 ± 5.9, which puts the eruption peak in 2010. 7±5. 9. The uncertainty on 
this postdicted date is not small, so the close coincidence between the postdiction and the 
real eruption date is fairly surprising. 

The error in the Schaefer (2005) and Schaefer et al. (2010a) predictions is simply the 
standard assumption that the light was being dominated by the usual thin a-disk (which 
yields E ^2). With the realization that the T Pyx light must be nonthermal and optically 
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thin, we can get E — 1 and an expected eruption around the year 2011. 

4. The Orbital Period Before Eruption 

In the classic review of RNe, Webbink et al. (1987) concluded that T Pyx must have a 
relatively short orbital period (like those of the typical classical novae) due to the lack of any 
red excess. Early photometric periods suggested were 0.069 days (Szkody & Feinswog 1988) 
and 0.1433 days (Barrera & Vogt 1989). With a much larger set of time series photometry, 
Schaefer (1990) found a stable and highly significant modulation but had severe daily alias 
problems, for which the highest peak in the Fourier transform was at 0.099 days, while 
another lower peak was at 0.076 days. Schaefer et al. (1992) collected an even larger 
set of time scries from 1966 to 1990 and determined the period to be 0.07616±0. 00017 
days. Patterson et al. (1998) made long runs from 1996 to 1997 to confirm this orbital 
period, discovered a steady period change, and discovered a complicated structure of weaker 
periodicities that vary in amplitude. Uthas et al. (2011) demonstrated that the stable 0.076 
day photometric period equals the orbital period as measured with a radial velocity curve. 
They also used many minimum timings from 1996-2009 (from the records of the Center for 
Backyard Astronomy) to give a high precision ephemeris that includes a parabolic term. This 
gives the minima times as 

HJD = 2451651.65255 + 0.076227249Ar + 2.546 x 10~^^A^^ (6) 

where N is the number of orbital cycles from the epoch. 

A highly accurate pre-eruption orbital period, together with a post-eruption period, will 
give the change in orbital period across the 2011 eruption. This period change will then give 
the ejected mass with high accuracy and no need for uncertain quantities (like the distance) 
or dubious assumptions (like the shell filling factors). For the case of T Pyx (but not for the 
other RNe), the dynamical friction of the companion star within the nova envelope will also 
cause a small period change (Livio 1991; Schaefer et al. 2010a) and this effect must also be 
taken into account. The ejected mass is critical for knowing whether the RN white dwarf is 
gaining or losing mass over the eruption cycle, and hence whether the RN will soon become 
a Type la supernova. Previously, only two other RNe (U Sco and CI Aql) had a measure of 
the period change across the eruption (Schaefer 2011). The T Pyx eruption of 2011 provides 
an opportunity to measure the period change for a third RN. 

As the period of T Pyx changes, it is vital to determine the pre-eruption orbital period 
for the time immediately preceding the eruption. To good accuracy. Equation 6 should work 
well to give the pre-eruption orbital period. With this, for an epoch of N=52662, the pre- 
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eruption orbital period is 0.07622993±0. 00000017 days. However, this period is based on 
minima times only until 2009. However, we also have the time series in February 2011, up 
to 48 days before the eruption. The possible existence of small period changes just before 
the eruption is critical where we need part-per-million accuracy. In this section, we will test 
the ephemeris of Uthas et al. (2011) and extend it from 2008 to February 2011 so as to get 
the period immediately before the eruption. 

With the time series photometry in Table 1, we have performed chi-square fits of the 
folded light curve to a template. The folded light curve is folded around the period from the 
parabolic cphcmcris of Equation 6. (Wc arc dealing with time scries of such short length that 
any approximate period would have been adequate.) The template we used was a flat light 
curve from phases 0.25-0.75 with a linear and symmetric decline to a minimum at zero phase 
of 0.15 mag fainter (see Figure 10 of Patterson et al. 1998). This procedure yields a formal 
uncertainty on the time of minimum, but this only accounts for the measurement errors. 
The uncertainty is actually dominated by the orbit to orbit variations in the light curve 
(see Figure 54 of Schaefer 2010 and Figure 8 of Patterson et al. 1998) caused by the usual 
flickering plus other unstable periodicities. When one orbit is considered, the light curve 
occasionally does not have an obvious minimum. However, when many orbits are considered 
together, the combined light curve always looks similar to the template. So the accuracy 
of the times of minima should scale as the inverse square root of the total duration of data 
in the time series in terms of the orbital period. For a group of similar observations over a 
relatively small time range, the scale of this uncertainty can be set by looking at the scatter 
in the 0-C curve about some smooth best fit line. So for example, the uncertainty in the 
minimum times reported by Patterson et al. (1998) can be calculated as 0.0033(A/0.076)~°'^ 
days, with A being the total duration of the observation (not necessarily contiguous) in days. 
With this. Table 4 presents the available times of minimum (in heliocentric Julian days) with 
their uncertainties. Included are six times of minimum from Patterson et al. (1998). 

These minima times can be compared to a linear ephemeris by means of the 0-C dia- 
gram. For the linear predictions, we merely take Equation 10 while zeroing out the quadratic 
term. For each minimum time, we calculate the number of cycles from the epoch (to the 
nearest integer), use equation 6 (without the quadratic term in this case) to calculate the 
model time, then difference this time from the observed time to form the usual 0-C curve. 
The 0-C values for the linear ephemeris are presented in Figure 2 as a function of the year. 
Uthas et al. (2011) reports that they have a large number of additional minimum times 
from 1998-2009, and these can well fill the gap in the 0-C curve. We see a nicely shaped 
parabola, which fits well with Equation 6. 

We have also sought deviations from the parabohc ephemeris. For this, we constructed 
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another 0-C curve where the ephemeris was Equation 6 (see the lower panel of Figure 2). 
If the Uthas et al. (2011) model is perfect, then all the observed minimum times should 
lie along the horizontal axis with O — C = to within the plotted error bars. Indeed, we 
see that the early time series from Schaefer et al. (1992), Patterson et al. (1998), and 
from early 2008 by Richards all fall close to the axis. This is to say that the parabolic 
ephemeris is a good match to the observed times from 1988-2008. The two times in 1986 
and 1987 are roughly 2-sigma and 1-sigma above the line, with a weak suggestion of deviation 
from Equation 10 before 1988, although we judge this suggestion to be too weak to accept. 
However, the two times from February 2011 are both closely consistent and with 3.7-sigma 
and 5. 1-sigma deviations from the parabolic ephemeris. These two measures indicate that T 
Pyx underwent a very small but significant change in period from 2008 to 2011. If the 0-C 
curve is linear from 2008 to 2011, then the orbital period at the time of the eruption was 
0.07622916±0.00000008 days. 

This result is encouraging, since the period has an accuracy of one part per million, 
which is in the right range to see the expected period change. However, this result is 
discouraging because it is 10 parts-per-million different from Equation 6. So it appears that 
the dominant uncertainty will likely be from very small changes in the orbital period, not 
from measurement errors. 

How long will it take to get a part-per-million post-eruption orbital period? With one 
full night of time series, the uncertainty in the minimum time should be roughly 0.0015 days. 
With two such observations separated by six years (with ~29,000 cycles), the uncertainty 
in the derived period will be close to one part-per-million. This duration can be improved 
by having more than just the two measures, with this improvement scaling as the square 
root of the number of minimum times. Over a three year interval, we could still get a 
part-per-million accuracy with eight full nights of photometry. 

5. The Rise to Peak 

T Pyx has a well-defined light curve with essentially continuous coverage and BVI colors 
starting at the the time of the discovery (Figure 3 and Table 2). The nova rose from V=13.0 
(discovery) to V=10.0 in a time of 8 hours, for an average rate of 9.0 mag day~^. This fast 
rise slows within the first day, rising to V=7.9, 2.0 days after discovery. Then, T Pyx fades 
by a third of a magnitude over the next day, followed by a very slow rise back to V=7.9 over 
the next five days. Soon after this, the rise increases to a minor peak of V=6.7 on JD2455685 
(19 days after discovery) and then up to the highest peak of Y—QA on JD2455693 (27 days 
after discovery). 



- 19 - 



We have measured the colors of the initial rise of a nova eruption (see Table 2). The 
B-V curve is plotted in Figure 4. Wc sec that the first measures (when T Pyx was around 
eleventh magnitude) have median B — V = —0.02, but this slowly changes to B — V = +0.23 
by two days after the start. At this time, just as the light curve is nearly fiat, the B — V color 
curve also goes nearly fiat. Just after discovery (around JD2455666.07), the median colors 
are B — V = —0.02 and V — I = +0.39. These colors are inconsistent with a blackbody, 
as for example an AO star with a temperature near 10,000K {B — V — 0.00) has the V-I 
color equal to 0.00, which is substantially different than that seen in T Pyx. Half a day 
after discovery (JD2455666.65), the median colors are S - F = +0.06, V-R^ +0.27, and 
V — I — +0.37. At the time when the fast initial rise breaks (JD2455668.0), the median 
colors are S - y = +0.23 and V - I ^ +0.58. 

The initial rise of a nova can be simply modeled by a uniformly expanding shell. The 
observed flux from the system will equal 

F^Fo+{vr/DfF,{T). (7) 

Here, Fq is some constant flux from the disk and companion star, v is the expansion velocity 
of the outer edge of the shell (corresponding to the effective photosphere), r is the time of 
expansion, D is the distance to the T Pyx system, T is the effective temperature of the 
photosphere, and Fjj(T) is the flux from a unit area of the photosphere. For the case of 
bolometric flux, F^(T) = aT* for a blackbody photosphere. Equation 7 can be converted 
into AB magnitudes as 

m = -48.60 - 2.5 logio[Fo + C{t - t^)], (8) 

with 

C={v/DfF,{T). (9) 

Here, r is expressed in the more practical form of the time {t) minus some fitted zero time 
for the start of the rise (to)- C is a combination of values that presumably remain nearly 
constant throughout the fast initial rise. Before to, the magnitude should be a constant 
corresponding to the flux Fo alone. We now have a simple model with three unknowns (Fo, 
C, and to) that should describe the initial rise. 

We have fitted our T Pyx light curve from the time of discovery until JD2455668.0 with 
a chi-square analysis. This includes 2151 magnitudes. To get a reduced chi-square of unity, 
we had to adopt a systematic uncertainty of 0.25 mag, which is a measure of systematic 
errors in the measurements plus deviations from the model. Our best fit model is shown in 
Figure 3. The constant flux corresponds to a magnitude of 14.63, to = 2455665.967, and 
C — 1.56 X 10~^^ erg s~^ cm~^ Hz~^. This model flts the initial fast rise well, and gives some 
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confidence that the simple model is close to the real case. The sudden deviation two days 
after the initial rise would then indicate some significant departure from the simple model 
in the sense that the expanding shell appears fainter than expected. We suggest that this 
deviation is caused by the outer portion of the ejecta expanding to become optically thin 
enough so that the photosphere can shrink significantly in size. 

Wc had originally hoped that this rise could be used to determine an accurate distance to 
T Pyx by the Baade-Wesselink method. (The distance to T Pyx is poorly known, for example 
with Schaefer (2010) giving 2500-4500 kpc.) Our light curve fit gives a value for C, the 
spectral line widths will give V, and we can use theory or colors to estimate Fy{T)). Then the 
distance can be directly derived from Equation 9. Unfortunately, the inputs are known only 
poorly. For the velocity, we do not know what position in a line profile (e.g., half- width- zero- 
intensity or half-width-half-maximum) corresponds to the velocity at the photosphere. For 
Fi,{T), the blackbody approximation is expected to be good for the photosphere, but the early 
colors are inconsistent with a blackbody. As such, the temperature cannot be determined 
with any useful confidence from observation of the colors. Also, no early spectrum exists 
so that we cannot get the temperature from spectral lines. From theory, the temperature 
varies by a factor of ~3, with a strong dependance on the white dwarf mass and other model 
details. With no way to get the temperature to any useful accuracy, we cannot get Fy{T) 
even if we know the emission mechanism. Nor can we turn this around to derive F^{T) from 
C, D, and V , because again the uncertainties of the inputs would be too large to be useful. 
In all, even though C is fairly well measured, we cannot get useful information on any of 
V, or F,{T). 



6. Accretion Rate 

We can use our data and analysis to address the question of the quiescent accretion 
rate. The value of M is key for understanding a variety of pecuharities in the T Pyx system. 
Perhaps the most critical question is reconciling the necessarily high accretion rate with the 
low X-ray luminosity. Here are four methods for determining M, some new, some old, and 
some failing. 

Method 1: The orbital period of T Pyx has been changing closely linearly with time 
(see the good parabolic shape of the 0-C curve in Figure 2). This large and steady period 
change can only arise from the mass transfer in the system. (Being below the period gap, 
T Pyx will not suffer from significant angular momentum losses from magnetic breaking, 
which can only make the period decrease.) The observed increase in the orbital period is 
expected when a low mass donor star transfers matter to the high mass white dwarf. The 
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period change will be proportional to M, so we can use the observed period change to deduce 
the accretion rate. The default scenario is conservative mass transfer, in which all the mass 
leaving the donor star arrives onto the white dwarf, and this material carries with it the 
specific angular momentum of the donor star. For conservative mass transfer, 

M = {-P/P){MwDM,omp)/[^{MwD - Mcomp)], (10) 

where P is the orbital period, P is the dimensionless time derivative of the period, 
is the mass of the white dwarf, and Mcomp is the mass of the companion star. The period 
derivative equals twice the quadratic coefficient divided by the linear coefficient (of the cycle 
count N) in Equation 6, so P = 6.68 x 10^^°. With Mwd = 1-3 M© (Schaefer et al. 2010a) 
and q = Mcomp/ Mwd = 0.20 ± 0.03 (Uthas et al. 2011), we have M = -3.5 x 10"^ Mq 
yr~^. (The negative value for M imphes that the companion star is losing mass, with the 
accretion rate often stated as the absolute value of this negative number.) 

How does the accretion rate change when the conservative mass transfer parameters are 
changed? For ordinary accretion disks, some amount of angular momentum must be lost 
so that material can move from the outer to inner parts of the disk, with this likely being 
carried by a few per cent of the material being lost through the L2 point. Warner (1978) 
derives the period change as a result of this 'liberal mass transfer' case, which when applied 
to the T Pyx case predicts M = —9.9 x 10"'^ Mq yr~^. Warner (1978) also points to another 
case of liberal mass transfer plus mass loss, where mass falling onto the white dwarf is then 
ejected from the system, a situation that would arise from nova ejecta or jets. For T Pyx 
in this case, M — —7.1 x 10~^ Mq yr~^. Patterson et al. (1998) point to another case for 
T Pyx, involving 'pure mass loss' where material is lost to the binary while carrying zero 
angular momentum. For T Pyx in this case, M — —2.5 x 10~^ Mq yr~^. We now have 
four cases with varying conditions, and the real situation is likely to be sufficiently complex 
(see the next Section) that none of these four cases is exactly correct. Nevertheless, these 
four cases demonstrate the range of possibilities and the range of required M values. From 
this, we see that the conservative mass transfer case yields the lowest accretion rate, with 
alternative possibilities perhaps raising the accretion rate by up to one order of magnitude. 
(The presence of any angular momentum loss and consequent period decrease would only 
force the accretion rate to larger values.) From the physics of the T Pyx situation, the best 
estimate is hkely the 'liberal mass transfer' case with M — —9.9 x 10"''' Mq yr~^, and we 
can say with some confidence that the conservative mass transfer case is the lower limit so 
the accretion rate is > 3.5 x 10~^ Mq yr"-*^. This result can be represented by an accretion 
rate of 10-6 0±o-5 Mq yr'^. 

Such an accretion rate appears to be very large, as it is near the limit for steady hydrogen 
burning (3 x 10~^ Mq yr~^ for a near Chandrasekhar mass white dwarf; Nomoto et al. 2007; 
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Shen & Bildsten 2007). This limit is for accretion uniformly over the entire white dwarf 
surface, so any channelling of the material onto a small area near the polar cap (as in the 
case of an intermediate polar like T Pyx) will result in a local accretion rate that will have 
steady hydrogen burning. For any plausible accretion rate and any plausible polar cap area, 
the local accretion rate will be higher than the limit for steady hydrogen burning, and this 
provides the luminosity that makes T Pyx a V1500 Cyg type star (Schaefer & Collazzi 2010). 

Method 2: The 2011 eruption ties the inter- erupt ion interval to 44 years. With this 
number, plus the approximate mass of the white dwarf, we can estimate the required trigger 
mass on the white dwarf and hence the accretion rate. A substantial problem is that the 
accumulated mass required for the trigger has a very strong dependence on the adopted 
conditions. It varies by over an order of magnitude over the acceptable range of white dwarf 
masses (Yaron et al. 2005), it varies by a factor of four over a plausible range of metaUicities, 
and it varies by a factor of five over a plausible range of ^He mass fraction (Shen & Bildsten 
(2007). In addition, the trigger mass varies by a factor of two with a change in M by a factor 
of ten over the plausible range. Nevertheless, for the estimated conditions, the T Pyx white 
dwarf must have accumulated lO-^-Oio.s (Yaron et al. 2005), 10-6-3±o.2 (Nomoto 
et al. 2007) or lO-5-5±o.3 ^^^^^ ^ Bildsten (2007) to be able to erupt. This range of 
theoretical models can be represented as io~^-^^°-^ Mq. With the 44 year interval, we get 
the average rate of deposition of mass onto the white dwarf is in the range iO~'''-^='='^-^ Mq 
yr~^. This number represents a lower limit on the M value (which is the mass lost by the 
companion star) because perhaps substantial fractions of the material in the disk could be 
ejected rather then accumulating on the white dwarf surface. Also the fraction of hydrogen 
in the accumulated material might be low due to steady hydrogen burning. Importantly, 
this highly uncertain estimate is based on specific theoretical models, so we must be careful 
to not use this estimated M to evaluate future models. 

Method 3: With our physical model for T Pyx (see Section 3), we can ask what accretion 
rate produces the observed brightness. That is, the system luminosity is tied to M, so we can 
work from the observed B=15.5 mag (valid for most of the time during the last pre-eruption 
interval) and the distance (3.5±1.0 kpc) and change the M value until we get agreement. For 
this, our model includes a disk extending from 50% to 90% of the Roche lobe and reprocessing 
of the white dwarf light. We need accretion rates of 10~^ °='=°-^ Mq yr~^ for the plausible range 
of distances and radii for the inner edge of the disk. SelveUi et al. (2003) report a similar 
calculation based on the ultraviolet luminosity, deriving accretion rates from 2.2 — 4.6 x 10~^ 
Mq yr~^. Selvclli et al. (2008) used a similar method to determine the accretion rates as 
1.1 ± 0.3 X 10~^ Mq yr~^, while they also report on various relations between M and the 
V-band absolute magnitude to come up with similar accretion rates. These calculations all 
assumed that the system light is dominated by ordinary blackbody emission from the disk 
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and companion star, whereas we have aheady seen that this cannot be. Therefore, these 
calculations are not relevant to T Pyx and the results can only be ignored. 

Method 4: T Pyx is at best only a faint X-ray source, and this can be used to place hmits 
on the accretion rate. Selvelli et al. (2008) report that their observations with XMM-Newton 
for 22.1 ksec in 2006 shows a faint source visible over the 0.2-8 keV energy range. The signal 

was too faint to allow for a reliable fit, but they demonstrate that a 200,000 K blackbody 
with 10^'' erg s^^ source at 3.5 kpc is rejected because the spectrum below 0.5 keV would be 
roughly 300 times brighter than observed. For the soft X-ray component, this would imply 
an X-ray luminosity of roughly lO'^'^"'"' erg s~^. Balman (2010) used the same XMM-Newton 
data set from 2006 (except she reports on effective exposure times of 30.6, 39.7, and 38.8 
ksec for the three detectors). She found that an X-ray image shows an extended profile, with 
the unstated fraction of light above a single point source certainly coming from the shell 
around T Pyx. A spectral fit gives a luminosity of 6.0 x 10^^ erg s~^ at a distance of 3.5 
kpc over 0.3-10.0 keV. If this X-ray luminosity is taken to be the light from the accretion 
onto the white dwarf {^GMwdM / Rwd), then the accretion rate is 8 x 10~^^ Mq yr~^. If 
the X-ray luminosity is taken to be from steady hydrogen burning of the freshly accreted 
material on the surface of the white dwarf (0.007X Mc^ if all the hydrogen is burned, where 
X is the usual hydrogen fraction), then the accretion rate is around 2 x 10^^^ Mq yr^^. 

With these estimates from four different methods, the obvious problem is the contra- 
diction between method 4 and all the other methods. The energetics of the 2011 eruption 
require a minimum mass to power the explosion and this was piled onto the white dwarf in 44 
years, and this places a firm lower limit on the average rate of accumulation onto the white 
dwarf (from method 2). The accretion rate in 2006 was not far from the average from 1967 
to 2011 (as measured by the B-band flux), so M must be greater than 2 x 10~^ Mq yr~"^, and 
likely is much greater than this. With this, the accretion energy onto the white dwarf must 
be > 2 X 10^^ erg s~^, which is 20 times larger than the value from Balman (2010). (Balman 
evaluated this discrepancy to be a factor of 10-100.) So in all cases, the X-ray luminosity is 
greatly lower than the minimum required by the accretion energy simply coming from the 
X-rays. With steady hydrogen burning on the white dwarf, this discrepancy becomes even 
larger. Therefore, we know that most of the accretion luminosity is not coming out above 
0.2 keV and therefore method 4 cannot be used to tell us the M. 

The question is now how to get the accretion (plus steady burning) luminosity to emit 
primarily below 0.2 keV. This can be done if the emission is from a blackbody with tem- 
perature <50,000 K. This fits well with the long-term ultraviolet spectrum from the lUE 
spacecraft (but not the optical or infrared spectrum) being well represented as a 34,000 K 
blackbody. Hence, the easy solution to the X-ray paucity observation is simply that the 
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accretion-plus-burning energy is coming out as a 34,000 K blackbody. 

What is the accretion rate in T Pyx? The first method (based on the P) is robust 
with 10~^°='=°-^ Mq yr~^. This purely dynamical method, derived only from Kepler's Law 
and the conservation of angular momentum, is based only on timings of high confidence and 
accuracy, while being independent of the usual problems (distance, extinction, and emission 
mechanism). The second method (based on the theoretical trigger mass) yields an accretion 
rate > IQ"''-^^''-^ Mq yr^-*^ and is also robust because it is based on simple energetics of the 
eruption, even though the exact value depends on the details of the model. The third method 
(reproducing B=15.5 mag from the disk) is unreliable because it assumes the usual thermal 
emission from disks and this cannot dominate the optical light for T Pjrx. The fourth method 
(based on the X-ray luminosity) is also not useful because it gives such low values for M 
that cannot allow T Pyx to be recurrent, so only some very small and unknown fraction of 
the accretion luminosity can be coming out above 0.2 keV. Only the first two methods are 
useful. We favor the first method because something close to the 'liberal mass transfer' case 
must be the situation for T Pyx, and we do not think that any other mechanism can produce 
such a good parabolic 0-C curve over two decades of time such that the orbital period is 
increasing. As a lower limit, the second method is in agreement with the first, although the 
limit is sufficiently far below the first method value that we expect that there must be some 
substantial inefficiency such as mass lost to the system from the accretion process or the 
partial burning of hydrogen on the surface of the white dwarf. So in all, we conclude that 
in recent decades the mass loss rate from the companion star has been 10~^ °='=°-^ Mq yr~^, 
while the white dwarf has accumulated material at the rate of > \{)-'^-^^^-^ Mq yr~^. 

With = 1, we can use the historical light curve to estimate the accretion rate at 
times outside the 1988-2008 interval. Before the 1890 eruption we have (Fb) = 7.59, while 
the 1988-2008 interval has (Fg) close to 1.47 mag (Table 3). This means that the early 
accretion rate was 7.59/1.47=5.2 times the value from the last paragraph. So around 1890, 
the companion star lost mass at the rate of io~^-^='='^-^ Mq yr~^. Another application is that 
the years 1978-2004 have an average B=15.5 mag, while from 2005-2011, T Pyx had an 
average B=15.7 mag. This corresponds to a difference of a factor of 16% in blue flux, and 
hence a difference in M by 16%. This provides an explanation why the P changes around 
2008 to a shghtly smaller value. That is, as the white dwarf luminosity secularly fades, the 
irradiation-driven mass loss from the companion decreases, and the orbital period change 
decreases. So we can work backwards from the historical light curve through the accretion 
rate to the orbital period as a function of time. With this, the fairly constant brightness of T 
Pyx from 2005-2011 makes for a fairly constant P over this interval. For a short time interval, 
like from 2008-2011, the curvature in the 0-C curve is negligible when simply extrapolating 
the orbital period from February 2011 to April 2011, therefore the pre-eruption orbital period 
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from Section 4 (0.07622916±0.00000008 days) is still good. 

7. The Nature of the Emission Prom Near the White Dwcirf 

T Pyx is a unique system, with its 1866 classical nova event that triggered the recurrent 
nova episodes, its very high accretion rate in a system with an orbital period of 1.83 hours, 
and its secular decline in brightness. There can be few stars in our Milky Way in this same 
state as T Pyx, but such rare systems have frequent eruptions to call attention to themselves 
from afar. T Pyx is also unique in the properties of the emission region near the white dwarf. 
The nature of this emission region has many possibilities and little theoretical understanding. 

We have put a variety of constraints on the central emission region. First, we know 
that the optical light is optically thin, because the observed inter-eruption intervals require 
Fb oc M which is impossible from a hot, optically thick source, and because the power 
law spectral energy distribution out to the infrared cannot come from any superposition 
of blackbodies within the Roche lobe. Second, we know that the accretion rate, i.e., the 
mass leaving the donor star, is roughly l0~^-°='=°-5 M© yr~-^ because the systematic and well- 
measured period increase from 1988-2008 cannot come from any cause other than the mass 
transfer, while the dynamical physics puts robust constraints on the accretion rate. Third, 
we know that the accumulation rate of material onto the white dwarf must be >iQ-'^-^^o-^ 
Mq yr~^, because this is required to provide the trigger mass and the energy for the 2011 
eruption. Fourth, we know that the region has a temperature that must be <50,000 K for an 
equivalent blackbody, because Balman (2010) used XMM-Newton to measure a limit on the 
X-ray luminosity of the central point source of 6.0 x 10^^ erg s~^ for photon energies >0.2 
keV. Fifth, the ultraviolet light is rising sharply towards the far ultraviolet as consistent with 
a 34,000 K blackbody and a luminosity of 2.9 x 10^^ erg s~^ (corrected to 3.5 kpc), as based 
on the /t/E spectral energy distribution (SelveUi et al. 2003). Sixth, the white dwarf has a 
magnetic field strong enough to disrupt at least the inner part of any disk and to channel 
the accretion onto a polar cap region, as shown by the unstable photometric periodicities 
near the orbital period reported by Patterson et al. (1998). 

Based on the strong precedent of other novae and cataclysmic variables, the default idea 
would be to have the flux coming from some sort of a boundary layer of an accretion disk 
or an accretion column plus additional light from an accretion disk with reprocessing. But 
given the uniqueness of T Pyx, there is little confidence in this stereotype. Indeed, we know 
that this conventional picture must be far from correct. The M value is around 10~^° Mq 
yr^^ being channeled onto a polar cap, so the default picture cannot explain the accretion 
region (Nomoto 1982). Also the power law spectral energy distribution in the optical and 
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infrared cannot be produced by any disk and the absence of a Rayleigh- Jeans tail requires 
non-thermal emission processes in the optical region. With the failure of the conventional 
picture, we must look to alternative possibilities. 

The mass faUing off the donor star has a lot of energy that will be radiated at some wave- 
length. The total luminosity from 10~^ ° M© yr"-*^ falling onto the white dwarf {GMwdM / Rwd) 
is 10^^'^ erg s~^. Not all of this energy needs to leave the system as radiation if a substan- 
tial fraction of the material falling towards the white dwarf is ultimately ejected. If only 
10^'^-^ M0 yr^-*^ makes it to the white dwarf, then the luminosity would be lO'^^''' erg s^^. 
The luminosity from complete hydrogen burning on the polar cap (0.007X Mc^) is lO"'*^'^ erg 
s~^ for the full accretion rate. The steady hydrogen burning cannot be complete because 
enough hydrogen must survive to make the nova eruption. In all, the accretion must provide 
a luminosity somewhere between 10^^-^ and 10^^-^ erg s~^. 

This accretion luminosity is easily identified with the flux seen by lUE from 120- 320 
nm with an apparent blackbody temperature of 34,000 K. This temperature explains why 
essentially no X-ray flux is visible with XMM. If this emission is blackbody at a temperature 
of 34,000 K, then the wavelength of peak intensity would be 85 nm, so only a fraction of the 
hght will appear in the lUE band. For a 34,000 K blackbody observed from 120-320 nm, the 
bolometric correction will be close to a factor of two. Applying this bolometric correction to 
the observed lUE luminosity gives 10"^^-^ erg s^^ as the luminosity of the central ultraviolet 
source. Thus, the observed lUE source can account for all of the accretion luminosity, if 
wc keep to the lower end of the range. Thus, there is little steady hydrogen burning and a 
substantial fraction of the material fails to reach the white dwarf. 

The luminosity and temperature points to the size of this ultraviolet-emitting region. 
For a 34,000 K blackbody equivalent with luminosity of 10^^-^ erg s"^, the surface area is that 
of a sphere with radius equal to 300,000 km. Such a volume is just smaller than the Roche 
lobe of the white dwarf. The coincidence of the volume required to produce the lUE source 
with the volume of the Roche lobe suggests that the binary orbit is somehow regulating the 
size of the emitting region. 

If the emission region nearly fills the Roche lobe, then the orbital motion of the com- 
panion star will stir the outer material causing substantial amounts to be ejected from the 
system. This could account for the difference between the mass lost by the donor star and 
that accumulating on the white dwarf. Such tidal effects will also provide a steady change 
to the orbital period, even though we have no clear idea as to the magnitude or direction of 
such period changes. 

It is unclear what fraction of this volume is optically thick. The inter-eruption intervals 
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give something close to Fb oc M, so that the region must be optically thin. But an optically 
thin region might not produce enough luminosity. The /t/E' spectrum only shows a continuum 
that follows Fi, oc u^-^^ (Selvelli et al. 2003), so this could either be a portion of an optically 
thick blackbody or an optically thin power law. Presumably, this volume will be denser 
towards the white dwarf, so there will be some radius outside of which it is optically thin, 
while inside some photosphere it will be optically thick. Further, the very hot accretion 
column near the white dwarf surface would likely be smothered by surrounding material 
that is optically thick in the inner regions, with this material suppressing the X-ray emission 
(Selvelli et al. 2008; Balman 2010) and also causing for the observed ionization state of the 
lines during quiescence (Gilmozzi & Selvelli 2007). 

The spectral energy distribution can be broken up into two components, the 34,000 
K blackbody that dominates in the ultraviolet plus the power law that dominates in 
the optical and infrared. The power law component is from optically thin emission, with a 
possible contribution from a circumbinary disk in the infrared. This suggests that T Pyx 
will be ejecting substantial amounts of material, perhaps into a circumbinary disk, which 
should be relatively cool and optically thin. All of the physical details (e.g., temperatures, 
densities, ionization states, and distribution in space) will depend critically on the nature of 
the hot inner region. Therefore, we can only assert that it is reasonable for the thin outer 
regions are optically thin and produce the power law spectrum throughout the optical and 
infrared bands. 

Our picture of the emission region (34,000 K, lO'^^ *^ erg s^^, the size of the Roche lobe) 
is without observational precedent among novae. Nevertheless, there is some theoretical an- 
ticipation that such an emission region could occur. In particular, when the accretion rate 
onto a white dwarf is between something like 10~^-^ to 10~^-^ Mq yr~^ (for a 1.3 M0 white 
dwarf with the accretion spread over the entire surface), a giant extended envelope will be 
formed around the white dwarf (Nomoto 1982). The T Pyx accretion rate (^iQ-^o±o-5 
yr~^) is easily within this range. Not all of the material necessarily makes it to the surface 
of the white dwarf (perhaps IQ^'^-^^''-^ M0 yr~^), but all this material will be concentrated 
by the magnetic field onto a small pole cap, so that the local accretion rate will exceed the 
rate to start making a large envelope. The physics of this situation will be very complex, 
and no theoretical study has yet been made of rapid accretion forming an envelope over only 
a fraction of the surface of a white dwarf. On top of this, there is also the complex effects 
of the nearby companion star stirring this whole envelope not far above the white dwarf 
surface. Any detailed theoretical understanding must have a full study involving magne- 
tohydrodynamics with a complex three-dimensional geometry involving nuclear burning on 
the surface, complex radiative transfer, and fast flows in many directions in a case where 
the key parameters (the accretion rate onto the white dwarf and the surface magnetic field 
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strength) have large uncertainties. We think that the situation is sufficiently complicated 
and uncertain so as to defy even simple theoretical expectations from having any useful con- 
fidence. So we are largely left with the observational picture for providing any useful idea of 
the actual physical state of the emission region. 

The observational picture of T Pyx shares some critical properties with one other cata- 
clysmic variable, V Sge (Patterson et al. 1998). The effective temperature is 27,000-40,000 
K and the optical-ultraviolet luminosity is iq^T-^^^^-^ erg s^^. The X-ray luminosity of V Sge 
is 3 X 10^^ erg s^^ from 0.1-3.5 keV with the preferred kT < 0.3 keV (Eracleous et al. 1991). 
V Sge has very strong blue colors and an ultraviolet excess, with B-V=-0.31 and U-B=-1.18 
after reddening corrections. V Sge has eclipses in its light curve whose breadth is nearly half 
the orbital period, hence indicating that the hot region is comparable in size to the Roche 
lobe of the white dwarf. These observed properties all require an emission region similar to 
that sen in T Pyx. Eclipse timings over a many year interval yield a parabolic 0-C curve 
with P/P = —5 X 10^ years (Patterson et al. 1998). V Sge has its period decreasing (whereas 
T Pyx has its period increasing), but this is easily explained by the high mass ratio g ~ 4 
(Herbig et al. 1965 give Mwj^ = 0.7 Mq and Mcomp = 2.8 Mq) so that the conservative 
mass transfer case implies an accretion rate of 10~^'^ Mq yr^^. This very high q (compared 
to all cataclysmic variables) readily explains why the accretion rate is so high, because any 
mass transfer only shrinks the Roche lobe of the donor star making for runaway accretion. 
With this, the cause for the high M in V Sge is completely different than in T Pyx. Another 
substantial difference between V Sge and T Pyx is that the V Sge white dwarf is far below 
the Chandrasekhar mass, and this explains why T Pyx shows frequent recurrent nova erup- 
tions while V Sge has had no observed nova events. Nevertheless, despite greatly different 
evolutionary paths, both T Pyx and V Sge apparently have hot sources with high luminosity 
that fill up the Roche lobe around a white dwarf caused by an accretion rate close to that 
required to initiate steady hydrogen burning. 



8. A Consistent Picture 

T Pyx has a fairly complicated history that has been pieced together from a wide variety 
of data. It might be useful to have one consistent picture explicitly collected in one place. 
Here, we present a largely chronological summary based primarily on this paper, Schaefer 
(2010), Schaefer & CoUazzi (2010), and Schaefer et al. (2010a). 

The basic system parameters of T Pyx are a white dwarf mass of 1.30±0.05 Mq, a 
companion star mass of 0.26±0.04 Mq, an orbital period of 1.83 hr, and a nearly pole-on 
inclination. The white dwarf has a sufficiently strong magnetic field so that the accreting 
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material will be channeled onto a polar cap, as with an intermediate polar. The distance is 
3.5±1.0 kpc, and the extinction has Eb-v — 0.25 ± 0.02 mag. 

For ~750,000 years before 1866, T Pyx was a nondescript short-period intermediate 
polar with an accretion rate of around 4 x 10~^^ Mq yr~^ and B — 18.5 mag. After this long 
and slow accumulation of mass on the white dwarf, a normal classical nova eruption occurred 
in the year 1866±5, reaching perhaps B — 5 mag. This eruption sent out a smooth shell with 
a mass of ~ 10^"^'^ Mq (larger than the accumulated mass on the white dwarf) at a velocity of 
500-715 km s^^. This eruption irradiated the nearby companion star (which is very close due 
to the short orbital period) so as to drive a high rate of mass being pushed through the LI 
point in the Roche lobe, with all this material being channeled by the white dwarf magnetic 
field into a small polar cap region. With some combination of the high accretion energy 
and steady hydrogen burning, the white dwarf remained as a hot high-luminosity source to 
continue the irradiation and keep the high M at a nearly self-sustaining rate. After the 1866 
eruption was long over, T Pyx stayed as a bright high-accretion system, much in contrast 
to its state a few years previously. With its newfound high accretion rate, the white dwarf 
accumulated material quickly, and with it being near the Chandrasekhar mass, the system 
started producing recurrent novae events. Each successive recurrent nova eruption sent out a 
fast (2000 km/s) shell of relatively low mass (~ 2 x 10~® Mq), with these shells ramming into 
the slower 1866 shell causing thousands of knots to form by Rayleigh- Taylor instabilities. 
The feedback loop of irradiation and accretion is not stable, so the accretion rate must be 
dechning, with the typical time scale for such cases being a century or so. With this, the 
accretion rate must dechne, the brightness of the system must decline (as observed to fade 
from B=13.8 in 1890 to B=15.7 in 2011), and the overall recurrent nova phase (starting in 
1866) must have a brief duration. Another consequence of the falling accretion rate is that 
the duration of the inter-eruption intervals must increase. 

Three requirements are needed for this history to apply to T Pyx: a white dwarf near 
the Chandrasekhar mass, a short orbital period, and a significant magnetic field on the white 
dwarf. The combination of all three is rare, so there will be few such stars in our Milky Way 
and the duration of their recurrent nova phase will be short, perhaps one-to-a-few centuries 
out of millions of years. Nevertheless, their frequent and bright eruptions will call attention 
to these systems from far across the galaxy. 

Although T Pyx is a system with material falling at a very high rate onto a white dwarf 
near the Chandrasekhar mass, it is not a progenitor of a Type la supernova. If the white 
dwarf is increasing in mass during the recurrent nova phase, stars like T Pyx cannot become 
the observed Type la supernovae because the former are by far too rare in our galaxy (and 
in the Andromeda Galaxy and the Large Magellanic Cloud) to account for the supernova 
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rate. Also the duration of the recurrent nova phase is such a small fraction of the system's 
lifetime that this would be much too low to account for the supernovae. Whether the mass of 
the white dwarf is increasing or not over the recurrent nova cycle is inconsequential because 
the dynamics are greatly dominated by the classical nova eruption and these always eject 
more material than was accreted, so the white dwarf is actually losing mass. In general, all 
recurrent novae with short orbital periods (only T Pyx and IM Nor are known as such) are 
not the supernova progenitors because they have no mechanism to sustain the high accretion 
rate. The long orbital period recurrent novae (e.g., U Sco and RS Oph) can still be supernova 
progenitors because they have a mechanism (the evolution of the companion star off the main 
sequence) to sustain a steady high accretion rate onto the white dwarf. 

The irradiation of the low mass companion star causes a high accretion rate, which is 
measured to be io~^ °='=°-5 Mq yr~^ as based on the fast and steady increase in the orbital 
period since the middle 1980's. With all this material being channelled onto a small polar 
cap, the local accretion rate will be so large that an extended envelope will form around 
the accretion column. The physics of the situation is very complex and has not yet been 
explored, but the observed result is a hot envelope that fills up the entire Roche lobe of 
the white dwarf. This hot envelope is seen in the ultraviolet by lUE as having a spectrum 
like a 34,000 K blackbody and a luminosity of lO'^^-^ erg s~^. The power law in the 
optical and infrared spectrum comes from the optically thin outer parts of the central region 
plus from the diffuse material outside the Roche lobe that has been ejected by the close-in 
companion star. With the optical light being from an optically thin region, the observed flux 
will be proportional to the mass accretion rate. As such, the 22 year inter-eruption interval 
before 1967 (with (B) = 14.70 mag) has a flux roughly a factor of two brighter than the 
inter-eruption interval after 1967 (with (B) = 15.59 mag), so its duration should be roughly 
2x22 years and the next eruption around the year 2011. 

The eruption after 1967 was widely expected starting in the mid-1980's. No eruption 
could have been missed (say, due to solar gaps) because the recurrent nova events have a long 
duration (roughly 230 days, with ts = 62 days) and a bright peak (B=6.7 mag) such that 
the nightly observations throughout the entire observing season would certainly have caught 
any eruption. Nevertheless, many observers in the professional and amateur communities 
persevered by taking many time series and nightly monitoring for decades on end. Prom 
1986-2011, the orbital period had a steady increase at P/P — -1-313,000 years, although a 
small deviation appears from 2008-2011 as expected for the small observed drop in accretion 
luminosity over this time. In the several years before 2011, T Pyx averaged V=15.5 mag, 
and was never seen brighter than V=15.0. 

It came as a surprise to Stubbings when T Pyx appeared at V=14.5 on UT 2011 April 
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5.51, and he realized that this might presage an eruption. Instead, he saw a unique rise and 
fall in brightness with a peak at V=14.4 over several days. This pre-eruption rise started just 
11 days before the eruption. Only four anticipatory rise/dips have ever been seen (T Pyx, , 
with all being greatly different from each other, and all with zero theoretical understanding. 

As part of his long series of monitoring of recurrent novae, Linnolt saw T Pyx at V=14.5 
on UT 2011 April 13.36 and realized that this was unusually bright. On the next night at UT 
2011 April 14.29, Linnolt saw T Pyx at V=13.0 and immediately recognized the eruption. 
Within 15 minutes the world had been alerted, Plummer and Kerr soon confirmed the 
eruption, and Nelson, Carstcns, and Streamer started time scries in BVl. For the first two 
days, the magnitude during this initial rise is closely represented as that of a uniformly 
expanding shell. The steady and frequent monitoring, principally by amateur observers, fast 
recognition and notification, and speedy response around the world all make for T Pyx being 
the first and only nova eruption where the light curve was followed throughout the entire 
initial rise. 
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Table 1. Pre-eruption Rises and Dip 



Nova 


Nova Class» 


^peak 




Vrtse 


AT (days)'= ^ 


Rise Morphology^ 


T Pyx 


P(62), RN, IP, short-Po^;,, shell 


6.4 


15.5 


14.4 


11 


'Bump' 


V533 Her 


S(43), CN, IP, shell 


3.0 


15.0 


13.3 


550 


Exponential rise 


V1500 Cyg 


S(4), neon CN, Polar, shell 


1.9 


20. 5^^ 


13.5 


>23 


Fast and bright rise 


T CrB 


8(6), RN, Red Giant donor 


2.5 


9.8 


10.8 


380 


Dip, variable depth and color 



''Data from Strope, Schaefer, & Henden (2010). 

''Class P has a plateau in the light curve. Light curve class S has a smooth monotonic decline. The parentheses contains 
the ts decline rate. CN indicates a classical nova with no known recurrences. IP indicates an intermediate polar system with 
a highly magnetized white dwarf star. 

'^Data from CoUazzi et al. (2009) and this paper. 

''at is the time from the beginning of the pre-eruption rise (or dip) to the start of the eruption. 



Table 2. T Pyx in quiescence (1890-2011) and the initial rise of tiie 2011 eruption 



Julian Date Band Magnitude Source 



2411491 


B 


13.80 


± 


0.15 


HCO 


plates 


2411493.6 


B 


13.80 


± 


0.30 


HCO 


plates 


2412211 


B 


14.10 


± 


0.15 


HCO 


plates 


2412236 


B 


14.30 


± 


0.15 


HCO 


plates 


2412936 


B 


14.30 


± 


0.15 


HCO 


plates 


2413570 


B 


14.30 


± 


0.15 


HCO 


plates 


2413891 


B 


14.30 


± 


0.15 


HCO 


plates 


2414424 


B 


14.50 


± 


0.15 


HCO 


plates 


2414454 


B 


14.70 


± 


0.15 


HCO 


plates 


2414766 


B 


14.70 


± 


0.15 


HCO 


plates 



(This table is available in its entirety in a 
ma<;hine-readable form in the online journal. 
A portion is shown here for guidance regarding 
its form and content.) 



Table 3. T Pyx Inter-eruption Brightness and Flux 



Interval AT (yr) Nights {B) (mag) (Fb) (Fb)AT (F^")AT 



Pre- 1890 




2 


13.80 


± 


0.13 


7.59 


± 


0.99 














1890-1902 


11.86 


12 


14.38 


± 


0.05 


4.52 


± 


0.20 


54 


± 


2 


65 


± 


3 


1902-1920 


17.92 


7 


14.74 




0.05 


3.20 


± 


0.15 


57 


± 


3 


67 


± 


4 


1920-1944 


24.62 


68 


14.88 


± 


0.03 


2.88 


± 


0.08 


71 


± 


2 


82 


± 


2 


1944-1967 


22.14 


28 


14.70 


± 


0.04 


3.35 


± 


0.11 


74 


± 


2 


87 


± 


3 


1967-2011 


44.33 


336 


15.59 




0.01 


1.47 


± 


0.011 


65 


± 


1 


69 


± 


1 
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Table 4. T Pyx Minimum Times 



Year Observer A (days) HJD minimum N Linear O-C (days) Quadratic O-C (days) 



1986. 


,022 


Warner'' 


0, 


,226 


2446439, 


,5121 


± 


0. 


.0035 


-68378 


0, 


.1267 


0.0073 


1987. 


.943 


Szkody'' 


0, 


,091 


2447141, 


,0770 


± 


0. 


.0056 


-59174 


0. 


.0960 


0.0065 


1988. 


.347 


Schaefer* 


1, 


,038 


2447288, 


,5651 


± 


0. 


.0017 


-57239 


0. 


.0844 


0.0006 


1988. 


,840 


Schaefer*' 


0, 


,154 


2447468, 


,9843 


± 


0. 


.0043 


-54872 


0. 


.0737 


-0.0033 


1989. 


,053 


Bond" 


0. 


,183 


2447546. 


,6591 




0, 


.0039 


-53853 


0, 


.0729 


-0.0012 


1989. 


,537 


Schacfcr'' 


0. 


,196 


2447723. 


,4983 


± 


0. 


.0038 


-51533 


0. 


.0649 


-0.0030 


1990. 


,342 


Buckley" 


0. 


,459 


2448017. 


,3506 


± 


0. 


.0025 


-47678 


0. 


.0611 


0.0030 


1990. 


,348 


Walker" 


0. 


,193 


2448019. 


,6358 


± 


0. 


.0038 


-47648 


0. 


.0595 


0.0014 


1996. 


,112 


Kempb 


3. 


,667 


2450124. 


,8300 


± 


0. 


.0005 


-20030 


0. 


.0095 


-0.0010 


1996. 


,220 


Kemp'' 


0. 


,892 


2450164. 


,3930 


± 


0, 


,0009 


-19511 


0, 


.0106 


0.0006 


1996. 


,289 


Kemp'' 


0. 


,328 


2450189. 


,5480 


± 


0, 


,0015 


-19181 


0, 


.0106 


0.0009 


1996. 


,352 


Kemp'' 


0. 


,386 


2450212. 


,4910 


± 


0, 


,0014 


-18880 


0, 


.0092 


-0.0002 


1997. 


,031 


Kemp'' 


1. 


,784 


2450460. 


,6090 


± 


0. 


,0007 


-15625 


0. 


.0075 


0.0010 


1997. 


,272 


Kempl" 


0. 


,686 


2450548. 


,4970 


± 


0. 


.0011 


-14472 


0. 


.0055 


-0.0001 


2008. 


,000 


Richards'' 


0. 


,160 


2454467. 


,1391 


± 


0. 


.0015 


36935 


0. 


.0334 


-0.0016 


2008. 


,019 


Richards'' 


0. 


,220 


2454474. 


,0788 


± 


0. 


.0013 


37026 


0. 


.0364 


0.0012 


2008. 


,058 


Richards'' 


0. 


,292 


2454488. 


,1785 


± 


0, 


.0011 


37211 


0, 


.0341 


-0.0015 


2011. 


,099 


Myers" 


0. 


118 


2455598. 


,7619 


± 





,0016 


51780 





,0627 


-0.0059 


2011 


il8 


Nelson' 


0, 


,297 


21.^)."')()0(i, 


, 00:59 


± 


0, 


,0011 


ol87o 


0, 


Ami. 


-(),()()."')7 



"Schaefer et al. (1992) 
''Patterson et al. (1998) 
"AAVSO data base 
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Fig. 1. — Pre-eruption rise. T Pyx never appeared brighter than V=15.0 for thousands 
of days before the eruption, and then it suddenly brightened to V=14.4 mag just 11 days 
before the eruption. (Further magnitudes from the All Sky Automated Survey are not shown 
here due to the apparent offset, yet they indicate that the duration is roughly 4 days for the 
pre-eruption rise.) The close time coincidence and the uniqueness of the bump demonstrate 
some causal connection between the bump and the eruption. The decline before the very 
fast rise demonstrates that this significant bump in the light curve is not some early part of 
the thermonuclear runaway. There is no understanding of the cause of the bump or why it 
anticipates the eruption. Explaining the pre-eruption phenomena is now a premier challenge 
for nova theory. 
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Fig. 2. — T Pyx 0-C curve. T Pyx has a stable periodic modulation in its light curve, and 
we have minimum times from 1986 until February 2011 (just 48 days before the eruption). 
The upper panel shows the 0-C curve for a linear ephemeris (Equation 6 with the quadratic 
term set to zero), in which we see a closely parabolic shape. This indicates that T Pyx has a 
steady change in orbital period, just as expected for steady conservative mass transfer at a 
rate of 3.5 x 10"^ M© yr~^, although this can only be a lower limit. The lower panel shows 
the 0-C curve for the predicted times coming from the quadratic ephemeris of Uthas et al. 
(2011). This steady period change well describes the minimum times from 1988 to 2008. But 
the two minimum times from February 2011 are consistent and highly significantly below 
the predictions of Equation 10, so apparently T Pyx underwent some sort of a small period 
change soon after 2008, such that its orbital period just before the 2011 eruption is slightly 
different (7.5 parts per million) from the period given by Equation 6. Such a change in P is 
expected due to the observed 16% decrease in system brightness (and hence a 16% decrease 
in accretion rate) from 2005-2011. 
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Fig. 3. — Initial rise. T Pyx now has the only nova light curve where the brightness 
was followed over the entire initial rise. For the first time, we see that the initial rise starts 
suddenly and is very fast. The rise is as fast as 9 magnitudes per day. A uniformly expanding 
shell is expected to rise fast in brightness followed by a substantial slowing (in magnitude 
units) after the first day, with the best fit model displayed by the curve in the figure. We 
see a sudden small fading and break from the model about two days after the start, with 
this perhaps indicating the recession of the photosphere. 
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Fig. 4. — B-V color over initial rise. We also have the first and only measure of the 
colors of the initial rise of a nova. The color does change somewhat over the first two days 
of the eruption (covered by the model fit in Figure 3), suggesting that the photospheric 
temperature does change over this interval, with this being contrary to the simple model of 
a uniform expanding shell with a constant temperature. The initial colors {B — V = —0.02 
and V — I — +0.39) are not those of a blackbody. 



